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Abstract. We show that the dynamical Wigner functions for noninteracting fermions and 
bosons can have complex singularity structures with a number of new solutions accompanying 
the usual mass-shell dispersion relations. These new shell solutions are shown to encode the 
information of the quantum coherence between particles and antiparticles, left and right moving 
chiral states and/or between different flavour states. Analogously to the usual derivation of the 
Boltzmann equation, we impose this extended phase space structure on the full interacting 
theory. This extension of the quasiparticle approximation gives rise to a self-consistent equation 
of motion for a density matrix that combines the quantum mechanical coherence evolution with 
a well defined collision integral giving rise to decoherence. Several applications of the method 
are given, for example to the coherent particle production, electroweak baryogenesis and study 
of decoherence and thermalization. 



1. Introduction 

Many problems in modern particle physics and cosmology require setting up transport equations 
for relativistic quantum systems in out-of-equilibrium conditions, including baryogenesis [1, 2, 3, 
4, 5] , heavy ion collisions [6] and out-of equilibrium particle production [7] . It is straightforward 
to write a formal solution for the problem in the Schwinger-Keldysh method, and detailed studies 
have been done on e.g. thermalization of quantum systems [8]. It is much more difficult to find 
a simple enough an approximation scheme that is usable in practical applications. Here we will 
report a recent progress [9, 10, 11, 12, 13] towards a scheme, in which one can account at least 
for certain particular types of nonlocal coherence effects in the presence of hard collisions. 

The key observation is that in cases with particular symmetries (the homogeneous time 
dependent case and stationary, planar symmetric case) the free 2-point Wigner functions have 
new singular shell solutions (at A:o = for the homogeneous system and at A;^ = in the static 
planar case), that carry the information about the nonlocal quantum coherence. When this 
shell structure is fed into the dynamical equations for the full interacting theory, they reduce to 
a closed set of quantum transport equations for the phase space distribution functions on the 
singular shells, which describe the on-shell particle numbers (mass-shells) and to a measure of 
the coherence (the new Uq^z = shells). 

The problem that initially led us to develop the current formalism was the electroweak 
baryogenesis by the CP-violating quantum reflection mechanism. We therefore briefly discuss 
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Figure 1. Shown are schematically 
how various quantities change over 
the profile of a planar phase transi- 
tion front. The symmetric (massless) 
phase is to the right and the broken 
(massive) phase to the left, and the 
front is moving to the right with a ve- 
locity Vuj- The green line presents the 
variation of the Higgs field. The red 
solid line shows the baryon number 
violation rate Tsph and the red dashed 
line marks the constant expansion 
rate H. The solid yellow line shows 
the left chiral quark asymmetry. 
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EWBG and the quantum reflection problem in this context. First, the baryon number is broken 
at nonperturbative level in the electroweak theory by the electroweak anomaly which couples 
the baryon and lepton currents to the Chern-Simons current in the gauge sector. The anomalous 
B-number violation is vanishingly small in the zero temperature [14] , but it becomes very fast at 
high temperatures. Above the electroweak phase transition temperature Tc « 100 GeV, where 
the electroweak symmetry is restored, the B-violation rate [15] Fsph = (20 — 25)a^T^ ~ 10~^r^, 
is much faster rate per degree of freedom than is the expansion rate of the universe during the 
transition: H ~ 10~^^T. Second, the P- and CP-parities, whose breaking is necessary to give 
direction to the B-violation, are broken by weak interactions themselves and by the complex 
phases in the quark mass matrix, respectively. CP-violating phases are also abundant (and 
needed) in many extensions of the standard model. A third condition for any successful BG- 
mechanism is to have out-of-equilibrium conditions. This is much harder to establish in the 
electroweak setting, because typical interaction rates between particles. Tint ~ ck^T" ~ 10^'^T, 
are enormous in comparison with the expansion rate during the transition. The only way to find 
out-of-equilibrium conditions with a standard expansion history is that the transition is of first 
order, proceeding through nucleation of broken symmetry (with (f) ^ 0, where (j) is the vacuum 
expectation value of the Higgs field) bubbles inside the symmetric = phase. In this case a new 
time scale, the bubble passing time, emerges: At;, ~ Iw/vw, where Iw is the width of the phase 
transition front and Vyj is the front velocity. Typically one finds Atb ~ lO^/T ~ FrJ. Because 
particles get their masses through the Higgs mechanism: m ~ ycj), fermions change from being 
massless to massive excitations as they cross the phase transition front, and this sudden change 
is enough to push particle distributions momentarily out of equilibrium. Careful studies of the 
EWPT dynamics have shown that bubbles are created far apart and they grow to enormous 
dimensions before their coalescence completes the transition. Thus, in electroweak baryogenesis 
one is led to study the CP- violating, out-of-equilibrium interactions of fermions with expanding, 
stationary and planar phase transition walls. This setting is illustrated in Fig. 1. 

The precise quantitative problem is to compute the total left chiral quark asymmetry in the 
front of the transition wall, as this asymmetry acts as the seed that biases the B-violating rate 
to create the eventual baryon asymmetry. This problem has been satisfactorily solved for the 
case of a very smooth phase transition front. This semiclassical baryogenesis formalism was first 



developed by use of WKB-techniques [3], and later from a more fundamental CTP-formalism [4]. 

However, when the wall is not particularly thick, the quantum reflection mechanism may become 
relevant. This problem is much harder to solve in the presence of interactions because of the 
inherent nonlocality of the reflection mechanism. However, it can be addressed by use of our 
coherent quasiparticle approximation (cQPA) scheme, and solving this problem remains one of 
our immediate goals. However, there arc also several other interesting applications for the cQPA 
formalism. These include a number of homogeneous problems, such as particle production by 
time- varying (oscillating) background fields, e.g. during the preheating after inflation, or at the 
end of phase transitions. We have not yet completed implementing the cQPA scheme to the 
EWBG-program, but we have applied it to several homogeneous problems, some of which will 
be illustrated in the sections below. 

This proceeding is organised as follows. In section 2 we introduce the basic CTP-formalism 
and introduce the Kadanoff-Baym equations in the mixed representation. In section 3 we derive 
the cQPA shell structure including the case with flavour mixing. In section 4 we find expressions 
for physical quantities and solve some simple collisionless problems to illustrate the physical 
content of the coherence solutions. In section 5 we consider the interacting theory and derive 
the density matrix equation with collisions for the homogeneous problems. In section 6 we extend 
our results for the scalar fields, and finally section 7 contains our conclusions and outlook. 



2. CTP-formalism 

In this chapter we very briefly review the Schwingcr-Kcldysh formulation of a relativistic 
fermionic quantum field theory. The mathematical quantity of interest is the 2-point Wightmann 
function: 

iG<f^{u,v) = {Mv)^^)) = T^{P M^)M'>^)} ■ (1) 

By definition is an ensemble average over a quantum density operator p, which fully describes 
the properties of the system. In the Schwinger-Dyson formalism one replaces the problem of 
finding p by a problem of finding all correlators defined by p. Since is an "in-in"-correlator, 
a complex time contour has to be introduced, one convenient choice of which is shown in Fig. 2. 
On such a Closed Time Path (CTP) [16, 17], the basic object is a path ordered 2-point function 
(Dirac indices are suppressed): 

iGc{u,v) = {Tc['ilj{u)i^{v)]), (2) 

where Tc defines time ordering along the Keldysh contour C. The correlator Gc{x,y) can be 
shown to obey a contour Schwinger-Dyson equation 

Gc{u, v) = GO(n, v)+ [ d^i / dSs G^ci^, zi)^c{zu Z2)Gc{z2, v) , (3) 
Jc Jc 

where the self energy function Sc couples Gq to an infinite hierarchy of equations for higher order 
Green's functions. The precise form of Sc depends on the model Lagrangian and a truncation 
scheme for the hierarchy. (In the weak coupling hmit a natural truncation is to substitute 
all higher than 2-point functions by their perturbative expressions.) The complex time valued 
function Gc{x,y) can be shown to obey a contour Schwinger-Dyson equation The real time is 
regained at the expense of introducing a 2 x 2 matrix structure; after some manipulations one 
can rewrite Eq. (3) in the following matrix form: 

G'-i®G = a3(5 + S®c73G, (4) 

where 

^-\G> Gp J ^ ^-\^> S^, 
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Figure 2. Schwinger-Keldysh 
path in complex time. 



and (T3 is the usual Pauli matrix, and we defined a shorthand notation ® for the convolution 
integral: 

f^g = [ d''zf{u,z)g{z,v). (6) 



Here Gp and Gp are the chronological (Feynman) and anti-chronological (anti-Feynman) Green's 
functions, respectively, and G^ and G^ are the Wightmann distribution functions, the former 
of which is just the correlator in Eq. (1). Above have also left out the labels u and v where 
obvious; for example 6 = 6'^{u — v). Self energy functions can be calculated for example by using 
the 2PI effective action techniques (see eg. [18, 17, 5]). For example: 

where T2 is the sum of all two particle irreducible vacuum graphs of the theory. Introducing 
the usual retarded and advanced propagators G'"(n, v) = 9{vP — v^){G'^ + G^) and G°'{u, v) = 
—6{v^ — u^)[G^ + G^), and a similar decomposition for the self energy function S, equations 
(4) can be recast into a particularly compact from: 

(Go ^ - S^'''') G^'" = 5 (8) 
(Go ^ - S") G<'> = S<'>®G^ (9) 

Equations (8) and (9) are called the pole equations and the Kadanoff-Baym equations [19], 
respectively. In general, the former will fix the spectral properties of the theory, while the latter 
will give the dynamical evolution, including quantum transport effects we are interested in. The 
KB-equation for G^ can be further rewritten as 

(Go ^ - Sh) G< - S< Gh = ^ (S> G< - S< G>) , (10) 

where the hermitian function Gh = ^ (G** + G"'). The equation for G^ is not needed, because 
G> and G< are related by the spectral function, A = \ {G^ - G") = |(G> + G<), and A 
and Gh are determined by the pole equations (8). The spectral function is further subject to 
a condition 2A{t,u;t,v)j^ = S^{u — v), which is just the direct space version of the famous 
spectral sum-rule and follows directly from the canonical equal time anticommutation relation: 
{ijj{t,u),'4>\t,v)} = -i5^{u-v). 

2.1. Mixed representation 

The next step, crucial to our program, is to move to the mixed representation Fourier 
transforming with respect to the relative coordinate u — v. However, it is necessary to first 
define more precisely the model we are interested in. With the electroweak baryogenesis and 
particle production problems in mind, we study a model described by the Lagrangian: 



C = i'4)^'il) + ipLfnipR + ipRm*ipL + £int , 



(11) 



where m{x) = mii{x) + imi{x) is complex, possibly spacetime dependent mass and >Cint is the 
interaction part to be defined later. All nontrivial interactions responsible for coherence, will be 
produced by the mass term. From Eq. (11) it is easy to find that 

G^\u,v) = 5\u - v){i -m*iv)PL - m{v)PR) , (12) 

where Pl,r = ^(1 T 7^) are the chiral projector matrices. Inserting the expression (12) into 
Eq. (10) and performing the Fourier (or Wigner) transformation we can rewrite the KB-equations 
in the mixed representation: 

(^ + i -mo - zm57')G< - e-^0{S^^}{G<} - e-'^{T,<}{GH} = CcoU, (13) 

where the G- and S functions are functions of the average coordinate x = ^(u + f ) and the k 
which is the conjugate momentum to r = u — v. For example the Wigner function G^{k,x) is 
defined as: 

G<{k,x) = j d^re'^-''G<{x + r/2,x-r/2). (14) 
The collision term is given by 

({S>}{G<}-{S<}{G>}) 

The C*-operator is the following generalization of the Poisson brackets: 

0{/}{5} = \ [dxf ■ dkg - dkf ■ dxg] 
and the mass operators rho and 7715 are defined as: 

mo,5G<{k,x) = mRj{x)e-'29^< G<{k,x) . (17) 

The KB-equation (13) together with the pole-equations (whose mixed representation forms we 
omit here) and the sum-rule form a complete set of equations for solving G< , A and Gh exactly, 
after the interactions, a scheme to compute S and the mass profiles are specified. In practice 
these equations are way too hard to be solved in their full generality, and several approximations 
are needed to find a solvable set of equations. The novelty of this work is to show that a tractable 
approximation scheme that is general enough to treat quantum coherence simultaneously with 
interactions docs exist. We shall now proceed to build our scheme by first constructing the full 
phase space structure of the free 2-point functions. 
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3. Free theory and the shell structure 

In the coUisionless case the last three terms in equation (13) are absent. However, the mass 
operators still contain infinite number of momentum derivatives acting on G'^. These gradients 
are controlled by space-time derivatives of equally high order acting on the mass function. So, 
in the regions where the mass is slowly varying we can truncate these expansions to the lowest 
order. Note that the spatial variations in G^ are not restricted by this approximation. In this 
limit we find 

{^+'^p^-mR-imn^)G< = 0. (18) 

Even this equation is very difficult to solve with general boundary conditions. We shall therefore 
restrict us to two special geometries relevant for the baryogenesis and particle productions 
problems: the spatially homogeneous and stationary planar symmetric cases. Consider first the 
homogeneous case. Observing that helicity is now a good quantum number, we can use the 
decomposition 

iG<^^^g<®\{l + hk-d) (19) 



Figure 3. Dispersion relations in the case of a scalar mass function. The dark filled (red) blobs 
show the mass-shell contributions and the light (yellow) blobs show the corresponding coherence 
solution at the new ko = shell in the homogeneous case (left) and at the kz = shell in the 
planar symmetric case (right). 



in the chiral Weyl representation. When this decomposition is inserted to Eq. (18), it eventually 
breaks into hermitian and antihermitian 2x2-matrix equations for g^: 

(H): 2ko9< = {H,g<}, (AH): idtg< = [H, g<] . (20) 

where H = — /i|k|(T3 + niRai — mja2 is immediately identified as the the local Hamiltonian. The 
algebraic H-equation constrains the form of possible solutions. Using the Bloch representation: 
9h ~ hido + 9i^^): it can be written as Babg^ = 0, which has nontrivial solutions only if 

det(S) = (A;g-k2- |m|2)A;g = 0. (21) 

Somewhat surprisingly, in addition to the usual mass-shell solutions there are new solutions 
living at shell A;o = (see Fig. (3)). The chiral structure of the mass-shell solution is 

9h,m-s = 2vrso fso {ko " ^|k|cT3 + mnai - mia2)6{k^ - \m^\) , (22) 

where sq = sgn(/co). This corresponds to the usual heHcity eigenstate that can also be found by 
the use of Dirac equation. The chiral structure of the ko = solutions is: 

5m-s = ^ {himnf^ - mift)as + \k\{f^ai + /2V2)) 6{ko) , (23) 

Here the density functions and some functions of time t and the three momentum 

k, whose evolution is determined by the dynamical equation. The mass-shell functions are 
related to the usual 1-particle densities of states and the new functions /{*2 describe the quantum 
coherence between particles and antiparticles. The complete solution to the algebraic equation 
is compactly written as 



9h='^^ 



(24) 



where = k^ -|- \m?\ and the matrices ph± and po can be read off from Eqs. (22) and (23). The 
full shell structure of the homogeneous model is shown in the left panel of the Fig. 3. 



A similar construction can be done for the stationary, planar symmetric case [9]. In that 
case the spin s orthogonal to the symmetry plane is conserved and eventually, in a particular 
frame [9], one finds equations very similar to Eqs. (20): 



isdzgf = 



Pg. 




(25) 



where P = koa^ — i{mjia2 + nijai) is the local momentum operator. The algebraic equation 
again gives rise to a singular shell structure; only this time the mass shell is accompanied by 
a singular momentum shell at k^ = 0, where k^ is the momentum perpendicular to the planar 
symmetry in the static rest frame. We do not write the explicit forms for these solutions here, 
but they can be found in ref. [9]. The complete shell structure for the correlator in the planar 
symmetric case in is shown in the right panel of the Fig. 3. 

The crucial issue leading to the new solutions in cither of the cases discussed above is that we 
did not assume translational invariance. Indeed, if one imposes translational invariance (in time 
in homogeneous and in space in the planar symmetric case), the dynamical equations in (20) 
and (25) become algebraic and it is then easy to show that they exclude the ko,z = solutions. 
Another important thing to observe is that the new ko,z = solutions do not appear in the 
spectral function despite the fact that in the collisionless limit the equation for the spectral 
function coincides with Eq. (18) for G^. The difference comes from the additional constraints 
from the spectral sum rule which are enough to kill the coherence solutions from A. Thus, as 
one might expect, coherence solutions exist only as correlations between particular 1-particle 
phase space states in the dynamical sector of the theory. 

3.1. Flavour mixing 

The discussion of this section can be straightforwardly extended to the case of several mixing 
fields. We shall consider here only the planar symmetric case. If we replace the simple 
complex mass m by a general complex mass matrix M, the analysis leading to equations (25) 
remains intact, apart from the fact that the local momentum operator becomes a matrix: 
P = koas~i{MHcr2 + MAHO'i), where M = Mh + iMah- The matrix structure makes extracting 
the spectral solutions more complicated. However, for a hermitian M the dispersion relation 
can still be expressed in a very simple form: 



where m^j refer to the mass eigenstates of M. This dispersion relation shows new interesting 
features. The diagonal entries obviously correspond to the simple DR's found above in the 
scalar mass case. The off-diagonal entries however correspond to the singular shells for the 
flavour mixing coherence. Interestingly, the latter DR's change continuously from one closely 
degenerate with mass shells and describing mixing between particles going mostly in the same 
direction to a shell close to /c^ = and describing the combined flavour mixing and reflection 
coherence between particles going on opposite directions. The complete phase space structure 
is shown in Fig. 4 for the case of a hermitian 2x2 mixing matrix. Clearly a very complicated 
set of mixing scenarios is possible with this dispersion relation complex, where one can have the 
two diagonal 1-particle densities coupling up to eight coherence functions. On the other hand, if 
the extrenous constraints are enough to restrict the mixing to, say, positive momenta only, but 
not enough to resolve the flavour mixing triplet (illustrated within the black circle in Fig. 4), 
one finds the usual flavour oscillation pattern widely studied in neutrino physics. 

The complex singular shell structure illustrated in this section is the foundation of, and gives 
the name to our approximation scheme: the coherent quasiparticle approximation or cQPA. 




(26) 




Figure 4. Shown is the complete set 
of dispersion relations corresponding 
to the case of 2 x 2 flavour mixing. 
The triplets of filled circles connected 
by dotted lines indicate states which 
at a given energy form the chiral 
mixing sub-complexes. The black 
filled circles indicate a possible flavour 
mixing complex when the chirality 
mixing is prohibited by extrenous 
constraints. 



4. Physical density matrix 

The dynamical equations in (20) and (25) look compelling. The time dependent equation (20) in 
particular has just the usual form of the evolution equation for a density matrix. However, they 
are actually meaningless because of the singular structure imposed by the algebraic constraints. 
That is, the 5^^-matrices are mere phase space densities that cannot be identified with physical 
density matrices. As was already suggested in the previous section, we will have to convolve 
them with some ko- and / or fc^-dependent weight functions that may quantify the prior extrenous 
information on the system [9] . Different types of such information functions were illustrated in 
the right panel of Fig. 3 by dashed lines. A physically observable density matrix can be formally 
written as 

pw{ko,kz,s;z) = 22 / ^iko,kz,s \\ kQ,k'^,s' ; z) g<,{k'^,ko;z) , (27) 

s' 

where W encodes the prior information. Let us make clear that this description is an effective 
one; in fact all restrictions should be imposed on the system directly by the collision integral. 
For example, our knowledge that a certain neutrino beam is moving in a given direction should 
be described by the time-dependent correlations imposed on the system by interactions with 
the apparatus that generates the beam. However, it is obvious that the effective description can 
often be more useful in practice. 

In many applications some parameters can be fixed by conservation laws while no prior 
information at all exists on some other parameters. For example, in our reflection problem 
the energy, the momentum parallel to the symmetry plane and the spin s are conserved, while 
we cannot assume any information on the perpendicular momentum of the state. The weight 
function appropriate for describing such a system then is W = lird^uj — ^0)5^,5'- Using this 
weight in Eq. (27) one finds that the density matrix is just the zeroth moment of the fundamental 
distribution with respect to kz'. pw = Psi^^, s) = {gf)- Note that in general there is no guarantee 
that the equation of motion for a weighted density matrix has a simple form, but in this particular 
case, because the momentum operator P depends only on ko, the collisionless equation for ps 
indeed is just 

isdzPs = Pps - PsP^ ■ (28) 

Unlike Eq. (25), the new evolution equation (28) is well defined because the matrix ps is a 
nonsingular object that depends only on the energy oj and spin s. A similar construction can 
be made in the context of a homogeneous problem, where one can fix the momentum k and 



helicity, and integrate over all energies. The resulting collisionless equation, now for a density 
matrix pw = /9^(|k|, h), is found to be the familiar one: 

idtPh = [H,ph], (29) 

where ph = Ph+ + Ph- + Po in terms of the matrices defined in equation (24). 



Moment connections 

Prom practical point of view the most important thing about the singular shell structure we 
have found is that it allows one to draw a one-to-one correspondence between the density 
matrix elements and the shell /j-functions. Indeed, hermitian p matrices contain four degrees of 
freedom, which is exactly the number of independent shell functions. Clearly, if the coherence 
shell solutions were ignored two degrees of freedom would be missing, making it impossible to 
describe coherence phenomena. Using the Bloch-representation p^ = ^((^o) + (ffi*)^*)' have 
the following relations between ((7^) (here the zeroth moments w.r.t. energy) and /±,i,2- 

(5o") = /+ + /- 

(5?) = ^(/^-/^) + /i' 

{9^2) = -^(/^-/^) + /2' 

(53') = -^?(/^-/^) + ^(^/Nl^/2), (30) 

where w = tuk. As we shall see later, these moment- /j function relations will allow one to 
rewrite the collision integral entirely in terms of the moment functions, leading to a solvable, 
closed system of equations for the density matrix. 



4-2. Physical quantities 

It is gratifying to see that the natural definition of the number density in terms of our mass-shell 
density functions coincides with some other nonequilibrium definitions found in the literature. 

Indeed, according to the usual Feynman-Stuckclbcrg interpretation the phase-space particle 
number density is identified as n = /+, while for (fermionic) antiparticles n = 1 — /_. With 
these identifications, using the inverse relations of Eq. (30) we find that: 

nx.h = ^{-hMgl)+mR{g1)-mi{9^))+\{g^) 

ny.h = ^{-h\k\{gl)+mR{g1)-mi{g!^))-l{g^) + l. (31) 

Setting Tr(p'*) = (^q) = 1 (physically this constraint corresponds to setting the chemical 
potential to zero) these expressions reduce to the ones obtained in ref. [20]^, where they were 
derived using the solutions to a Dirac equation and a Bogolyubov transformation to diagonalize 
the fermionic Hamiltonian. 

One can also compute the energy and the pressure densities from the diagonal elements of 
the energy moment tensor 9'^'^ = ^{^^'^d'^'^ — d'^^^'^'^) + p ^ v. For energy density one finds 

{n{x)) = {9'\x)) =Y^j —u^^ (nk, + n^u - 1) , (32) 



^ Different signs of terms involving h and mi are due to a different sign convention in our definition of the 
Hermitian Wightmann function. 



which, according to expectations is simply a sum of free particle and antiparticle contributions, 
plus the vacuum energy. For the pressure we find 

= (r (x)) = E / 3 ("k. + n^h - 1) - mnf^ + rujf^^ . (33) 

Now, in addition to normal free particle and antiparticle terms there is an explicit contribution 
from the coherence shell functions /i,2- Thus the pressure is different from the statistical one 
at the quantum level. However, the coherence functions /i,2 are oscillatory with microscopic 
time-scales Atosc ~ 1/w and the classical thermodynamical pressure arises when Eq. (33) is 
averaged out over time-scales exceeding the quantum scale Aiosc- 

4-3. Solutions: coUisionless examples 

We applied Eq. (28) to the Klein problem of reflection off a step potential wall in ref. [9] . Correct 
results were found for the transmission and reflection coefficients and tunneling factors when the 
coherence shell factors are included in Eqs. (30), while neglecting them gave only pure classical 
complete transmission and complete reflection solutions. Similar results were later found also for 
the bosonic Klein problem [10]. We also considered a reflection off a smooth mass wall, shown 
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Figure 5. Reflection from a potential 
wall due to spatially varying complex 
mass function. 



in Fig. 4. This is precisely the problem of interest for the quantum reflection mechanism in the 
electroweak baryogenesis. Again, the integrated moment equations are exact for this case, and 
the moment connections are needed only asymptotically to deflne the reflection coefficients. The 



Figure 6. Shown are the mass-shell 
functions f± corresponding to the left 
and right moving fluxes and functions 
/i,2) which encode the quantum 
coherence. We have taken s = 1 and 
g/|m_oo| = 0.088. 
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solutions for the mass-shell functions and /+ (here equivalent to the left and right moving 
currents, respectively) and the coherence functions /i,2 are shown in flgure 6. The coherence 
functions are oscillating in the massless phase and vanish in the massive phase. This is as 
expected because the boundary conditions (see Fig. 4) exclude right moving states at z — > — oo. 
Our results for the reflection and transmission currents agree with those found earlier with use 
of the Dirac equation [21]. The crucial improvement over these results is to include the effects 




of interactions. We have not yet completed this work for the planar symmetric case. However, 
parallel results for the homogeneous problem do exist and they will be presented next. 

5. Interacting fields 

With interactions included the full Schwinger-Keldysh equations become much more complicated 
of course. A tractable approximation scheme can still be obtained by combining the 
familiar quasiparticle approximation with our zeroth gradient order treatment of the constraint 
equations. In this scheme we can in particular drop [10] the G//-term in Eq. (13). Eventually 
we find: 

(1^ + i -mo - zmsT' - ^h)G< = Ceoii • (34) 

Here we are focusing only on the effects of the nontrivial mass function, so we drop the thermal 
corrections to the dispersion relation setting 'Eh = 0. After inserting the helicity decomposition 
Eq. (19), breaking the equation to the hermitian and antihermitian parts and, consistently 
with our cQPA scheme, dropping the collision integral in the constraint equation, which is also 
expanded to the lowest order in gradients^, we find: 

(H): 2kog^ = {H,9<} 

(AH): dt9^ = -i[H,g<]-Cj;, (35) 

where is the antisymmetric part of the helicity projected collision integral. In cQPA and for 
an uncorrelated (for example a thermal) self-energy it can be written as [13] 

iC^ = cos{ldld!){r^,g<} + zsmCAdf)[Tn,g<] - {r^,g<^} , (36) 

where Th = |(S^ + i;^) and we defined the hermitian helicity projected self energy function 
Phj^T,'^'-^ Ph = S^'^ ^{1 + hk ■ a) . Finally g^^ is the equilibrium limit of the distribution g'f^. 
We were careful to keep a subset of gradients to infinite order, because the projection onto ko 
shell of these expansions is not controlled by any small parameters (F or m') and need to be 
resummed. Of course equation (35) is equally badly defined as was our free theory evolution 
equation for g"^. It is also clear that when the integration is done, the integrated collision term 
will, thanks to the moment connections Eq. (30), reduce to and expression entirely expressible 
in terms of the density matrix components. When integration (over ko with a fiat weight) and 
resummation is done, one eventually finds [13]: 

dtPh = -iW, Ph] - {r/is, 5ph} - {Tha, Sph+ - Sph-} - ^ [Tha, [H, Ph\] ■ (37) 

Here brackets {, } denote an anticommutator, ph = ph+ + ph- + Po in terms of the matrices in 
Eq. (24) and Spi = pi — pi^eq, and the collision terms F/j(s = ^(F/i+ it Th-) denote symmetric 
and antisymmetric combinations of the mass-shell interaction strengths F^^ = Th^ko = ztojk). 
Equation (37) is the sought after generalization of the standard collisionless equation of motion 
(29) for the density matrix to the case with decohering collisions. 

We have applied the equation of motion (37) to study a system with a temporally evolving 
CP- violating mass function: 

\m{t)\ = ^(1 +tanh(t/r^)) , arg[m(t)] = ^Ae{l - tanh(t/r.^)) , (38) 

^ Note that these are just the standard assumptions that one makes when deriving the usual Boltzmann equation 
from the Kadanoff-Baym equations. See e.g. [22]. 




Figure 7. Shown are the total asymmetry A(?7,) = X]/i(^k/i. — "^k/i) (see Eq. (31)) and the 
L-chiral current density per unit k volume: ji]^{t) = YlhC^ihi^ + ^^)Ph] ~ 5) (vacuum part 
subtracted). Blue curves correspond to the noninteracting and black curves to the interacting 
cases. 



where Tw is the characteristic time of the transition. Collisions were modelled by an interaction 
with a thermally equilibrated fermion and scalar fields qr and 0: 

Ant = -y ^l4> Qr + h.c. (39) 
The lowest order contribution to the self energy from this interaction is: 

^<'>{k,x) = i\y\'l ^^P^G<^>{k',x)PLA<'>{k - k',x) . (40) 

Explicit analytical forms for S*^'^ can be found in ref. [10]. Numerical results for the total 
asymmetry and left chiral current density are shown in Fig. 7. For various parameters we used 
(in units of |k|) = 2, A9 = ^, T = 10, ruq = 10, rUfp = 5 and y = 0.3 (0 in the free case) and 
the initial conditions at t ^ —00 were set to to thermal equilibrium: fl^ — > l/(e^^'^'= + 1) and 
fi2 ^0. In the collisionless case no total asymmetry is created and the L-chiral asymmetry 
oscillates with a constant amplitude after the transition. In the interacting case a temporary 
nonzero total asymmetry is created and both the total asymmetry and the amplitude of the 
oscillating L-chiral asymmetry decay as a function of time. Such a simple picture might be used 
to model the baryon asymmetry generation during a second order (or cross-over) electroweak 
phase transition, for example in models where the expansion rate of the universe is highly 
accelerated during the EWPT. 

6. Scalar fields 

Our method can also be applied to the bosonic fields with similar results. In the case of the 
scalar fields we wish to study the correlator: iA'^{u,v) = {(p{v)(f){u)) = Tr{p (j){v)(j){u)}. The 
derivation of the KB-equations for proceeds analogously to the fermionic case [11]; the final 
result in the mixed representation is 

Aq iA< - e-'0{UH}{A<} - e-'0{U<}{AH} = Ccoii , (41) 

with 

Ccou = -ie-'^ ({r}{A<} - {n<}{^4) , (42) 



where we already used the connection = |(A^ + A<) and defined T = |(n> + n<). The 
self energy functions are again computable from the 2PI-generating function: Il'^'^{u,v) = 
—i5T2[A]/6A^''^{v, u). Of course the explicit forms of Aq as well as the interactions depend on 
the model. Similarly to the fermionic case, we are interested in effects caused by a space-time 
varying mass function, so we take 

C = l{d^(Pf - \m{xfct? + Ant , (43) 

where >Cint is some yet to be defined interaction term. The inverse free propagator corresponding 
to Eq. (43) in the mixed representation then is 

A-i = k'-\d'' + ik-d- m\-'^^^< , (44) 

where the 5^-derivative always acts on the mass term and the 5j^-derivative to A^. We again 
consider the free theory first, and also expand the Aq operator to the lowest nontrivial order in 
gradients. Because the operator is complex, we find two independent equations: 

(kl-k^ -m'^{t) - \d^yA<{k,t) = (45) 
kodtiA<{k,t) = 0. (46) 

Unlike with fermions, the nonintcracting scalar KB-equations can not be divided into algebraic 
constraints and dynamic equations because both equations (45-46) contain time derivatives. Yet 
a similar shell structure does emerge. Let us first assume that fco / 0. Then^ Eq. (46) requires 
that dtiA^ = at all times and so one must also have dfiA'^ = 0. Substituting this back to 
Eq. (45) gives 

{kl-k^-m\t))iA<_, = 0. (47) 
This equation has the spectral mass-shell solution parametrized by t: 

iA<_,{ko, Ht) = 27rsgn{ko)fs,^{\k\,t)S{kl - - m\t)) , (48) 

where sq = sgn(fco). However, the point feg = is special: if we first set ko = 0, then equation 
(46) is identically satisfied, but Eq. (45) becomes 5|A< = — 4a;^(t)A^. For a generic time- 
varying mass this implies a corresponding spectral solution for kg = 0: 

iA<_^{ko, |k|,t) = 27T U\k\,t)S{ko) , (49) 

where /c(|k|, t) is some real-valued function. The most complete solution for a given momentum 
|k| is the combination of the solutions (48) and (49): A< = A^_g -|- AQ_g. The situation is now 
seen to be qualitatively equivalent to the case with fermions and we interpret analogously that the 
new ko = 0-solution (49) describes the quantum coherence between particles and antiparticles. 

6.1. Moment connection 

Let us again suppose that we do not have any prior restrictions on the momentum variable k. It 
is then appropriate to use the n-th moments as the weighted density functions for the problem: 

Pn{\k\,t) = J ^kEiA<{ko,\k\,t). (50) 

^ Note that including collisions at this stage would destroy the singular phase space structure. Neglecting them 
here corresponds to the usual quasiparticle approximation. 
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Figure 8. Shown is the number density (thick soUd hne) and the coherence function /c(|k|) 
(dotted(blue) hne). For the driving mass function Eq. (53) we used the parameters = T, 
A = 1.5T and B = O.IT. For other parameters we used |k| = 0.6T, y = 1, = O.IT and 
oj^f, = O.IT, where temperature T sets the scale. At r = co^pt = the system was in the adiabatic 
vacuum, and the driving mass term smoothly set to constant m = 2.5 T after r > 47r. 



When one computes three lowest moments using the complete shell structure of one 
immediately finds: 

PI = Uf+ + f~) 

P2 = ^(/+-/-). (51) 

All higher moments can be trivially related to these three lowest ones. As was the case with 
fermions, these relations are only valid to the lowest order in mass gradients. However, when 
used as an ansatz in the complete set of equations (up to all order in gradients in the flow term) , 
they lead to a closure and a solvable set of equations for the moment functions pn- 

6.2. Evolution equations with collisions 

When the cQPA ansatz is inserted to the full KB-equation (41) (again dropping terms with 
IS.H and n^f), and the resulting equations are integrated over the energy, and the projections 
of the collision terms onto /cq = shell are resummed, one finds the following closed set of 
equations [23]: 

\d'fpQ + ujIpq - P2 = - 2t^Tm dtpo 

^tPi = -^^m{pi - pl,cci) 
dtP2-\dt{m^)pO = -^Tm{p2- P2,eq) ■ (52) 

We illustrate the solution of Eqs. (52) in the case where the scalar mass term is an oscillatory 
function coming from a coupling of (j) to an oscillating background field ip (the inflaton): 

m^{t) = \mQ + Acos{2uj^t) + iBs\n{2uj^t)\'^ , (53) 

Physically this situation could model the particle production after inflation (preheating). 
Interactions were modelled by decays and inverse decays caused by the interaction with a 
thermalized fermion background: £int = where is some generic thermalized background 

fermion field. As the inflaton oscillates around its minimum, (/)-fields are excited out from the 
temporally varying classical background and the number density increases in each oscillation 
cycle. Also the magnitude of the coherence function increases at each oscillation, although the 
decohering effect of collisions is clearly visible between the oscillation periods. After the driving 
mass terms is set to constant, a phase of smooth decoherence and thermalization is observed. 



7. Conclusions 

We have introduced a novel approximation for the KadanofF-Baym equations for the fermionic 
and bosonic 2-point functions. Our approach makes use of newly discovered extended singular 
phase space structure for these functions in the free theory limit in two simple geometries: the 
spatially homogeneous and stationary planar symmetric cases. When applied to the full KB- 
equations, this coherent quasiparticle approximaUon (cQPA) leads to equations of motion for 
a density matrix (or a set of moment functions for scalar fields), that can describe nonlocal 
coherent evolution of the state in the presence of hard collisions. In the homogeneous case the 
new shell solutions characterize the particle-antiparticle correlations and in the planar symmetric 
case correlations between reflecting states moving on opposite directions. Our original, and 
still principal motivation for developing this formalism was the desire to find a quantitative 
description for the electroweak baryogenesis proceeding through the CP-violating reflection 
mechanism. While we yet have to complete this program, we have already demonstrated the 
use of our formalism in a number of other interesting problems; for example the homogeneous 
time dependent problem of particle production in the early universe at the end of inflation or 
sonic phase transition. This formalism can be used to model homogeneous baryogenesis as well. 
Yet another interesting research direction is to use the cQPA to better understand the nature 
of the flavour coherence in the case of neutrino mixing and in particular in connection with the 
supernova neutrinos. 
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